Analogous Hawking Radiation from Astrophysical Black Hole Accretion 
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For the first time in literature, we provide an exact analytical scheme for calculating the analogous 
Hawking temperature and surface gravity for general relativistic accretion onto astrophysical black 
holes. Our calculation may bridge the gap between the theory of transonic astrophysical accretion 
and the theory of analogous Hawking radiation, and may allow one to compare the properties of the 
two types of horizons, electromagnetic and acoustic, in connection to the black hole physics. We 
show that the domination of the analogous Hawking radiation over the actual Hawking radiation may 
be a real astrophysical phenomena. We also discuss the possibilities of the emergence of analogous 
white holes around astrophysical black holes. Our calculation is general enough to accommodate 
accreting black holes with any mass. 
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By providing an analogy between the propagation of 
light in a curved geometry and the acoustic propagation 
in a steady, inviscid, transonic barotropic fluid, it has 
been demonstrated [l| that the representative fluid posses 
the signature of an acoustic black hole (BH hereafter) 
embedded by an acoustic horizon (AH hereafter) at its 
transonic point (where the fluid Mach number becomes 
unity). This AH emits thermal radiations known as anal- 
ogous Hawking radiation ( AHR hereafter) . The tempera- 
ture of such radiation is called Analogous Hawking tem- 
perature, or acoustic Hawking temperature (Unruh |l| 
was the first to propose and compute such temperature). 
Throughout our work, we use Tah to denote the analo- 
gous Hawking temperature, and Th to denote the the ac- 
tual Hawking temperature where Tjj = tic 3 /SirkbGMsH, 
Mbh being the mass of the black hole. We use the words 
'analogous', 'acoustic' and 'sonic' synonymously in de- 
scribing the horizons or black holes. Theory of AHR 
opens up the possibility to experimentally verify some 
basic features of BH physics by creating the sonic hori- 
zons in the laboratory. A number of works have been 
carried out to formulate the condensed matter or optical 
analogous of event horizons |2j • It is also expected that 
AHR may find important use in the field of Inflationary 
models, quantum gravity and sub-Planckian models of 
string theory Q ■ 

Since the publication of the seminal paper by Bondi in 
1952 the transonic behaviour of the accreting fluid 
onto compact astrophysical objects have been extensively 
studied in the astrophysics community (including the fast 
ever paper on curved acoustic geometry by Moncrief 5]), 
and the pioneering work by Unruh in 1981 0] initiated 
a substantial number of work in the theory of AHR with 
diverse fields of application stated above. However, it is 
quite astonishing that, any attempt is yet to be made to 
make a bridge between these two categories of research, 
the astrophysical BH accretion and the theory of AHR, 
by providing a self-consistent study of AHR for real as- 



trophysical fluid flows. Since both the theory of tran- 
sonic astrophysical accretion (TTAA hereafter) and the 
theory of AHR stems out from almost exactly the same 
physics, the propagation of transonic fluid with acoustic 
disturbances embedded into it, we feel the pressing need 
to study AHR for relativistic transonic accretion onto 
astrophysical black holes and to compute Tah for such 
accretion. 

We consider general relativistic (GR) spherically sym- 
metric accretion of inviscid polytropic fluid onto a 
Schwarzschild BH of mass Mbh ■ We first formulate and 
solve the conservation equations for such accretion and 
then demonstrate that such flow may become transonic 
depending on initial boundary conditions. We then ana- 
lytically calculate the exact location of the acoustic hori- 
zon Th |6| and calculate the quantities on AH which are 
required to compute the surface gravity and Tah for our 
model. We then calculate Tah in curved space time for 
all possible real physical accretion solutions. 
The importance of our work, as we believe, is not only 
lies in the fact that we, for the first time, completely 
analytically compute the exact value of Tah for GR as- 
trophysical fluid flow in curved space time |lflj . Our cal- 
culation rather bridges the long standing gap between 
two apparently different school of thoughts sharing in- 
trinsically the same physical origin for their work. And 
perhaps the most important point to note here that the 
accreting astrophysical BHs are the only real physical 
candidates for which both the horizons, actual (electro- 
magnetic) and analogous (sonic), may exist together si- 
multaneously. Hence our application of AHR to TTAA 
may provide a novel procedure to compare the proper- 
ties of these two kind of horizons in its most consistent 
way i.e., by studying and comparing the behaviour of the 
same flow close to these horizons. As an illustrative ex- 
ample of such procedure, we also calculate and compare 
Tah and Th for the same BH and for same set of initial 
boundary conditions, and show that for astrophysically 
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interesting regions of parameter space, Tah may well ex- 
ceed Th which indicates that sometimes AHR may be 
dominant over the actual Hawking radiation phenomena. 
Also we find the possibility that an actual BH event hori- 
zon may be embedded by an acoustic white hole. 
We take the Schwarzschild radius r g = 2GMbh/c 2 , and 
scale the radial coordinate r in units of r g , any velocity 
in units of c, and all other physical quantities are scaled 
accordingly. G = c = Mbh = 1 is set. Mass of the 
accreting fluid is assumed to be quite less compared to 
Mbh (which is reality for astrophysical BH accretion), 
so that the gravity field is controlled essentially by Mbh 
only. Accretion is 9 and </> symmetric and posses only 
radial inflow velocity. We concentrate on stationary so- 
lutions. Accretion is governed by the radial part of the 
GR time independent Euler and continuity equations in 
Schwarzschild metric. The conserved specific flow energy 
£ (the relativistic analogue of Bernoulli's constant) along 
each stream line reads £ = hut where h and u^ are the 
specific enthalpy and the four velocity, which can be re- 
casted in terms of radial three velocity u and polytropic 
sound speed a to obtain: 
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7 — (1 + a 2 ) J y 1 — u 



(1) 



Substitution of Uh and ah in eq. (1) for r = rh provides: 



where 



rlTx + r h T 2 +T 3 = Q 



2£ 2 (2-3 7 ) + 9( 7 - 1) 



(6) 



4(7-1) (£2-1) 



r 2 = 



£ 2 (3 7 - 2) - 27 (7- 1)' 
32 (£ 2 - 1)( 7 - l) 2 



, r 3 = 



27 



64 (£2-1) 
(7) 

Solution of eq. (6) provides the location of AH in terms 
of only two accretion parameters {£, 7}, which is our two 
parameter input set to study the flow, and hereafter will 
be denoted by Vi- We solve eq. (6) completely ana- 
lytically by employing the Cardano-Tartaglia-del Ferro 
technique. We define: 



Si = 



3r 2 - r? _ 91^2 - 27r 3 - 2T\ 



£2 = 
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The mass accretion rate M may be obtained by integrat- 
ing the continuity equation: 



M = An pur 2 



1 



r (1 — u 2 ) 



where p is the proper mass density, p = Kp 1 equation of 
state is used (use of which is common in astrophysics to 
describe relativistic accretion) to define another constant 



a as: 



3 = K 1-7 M = Air pur 
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a 2 (7 - 1) 



7- (l + a 2 ) 



(3) 



Simultaneous solution of eq. (1-3) provides the dynami- 
cal three velocity gradient at any radial distance r: 



du 
dr 



u(l- u 2 ) [a 2 (4r - 3) - l] 



2r (r - 1) 



The denominator of eq. (4) becomes zero at r = 1 (actual 
horizon) and when u = a (acoustic horizon). To main- 
tain the physical smoothness of the flow for any r > 1, 
the numerator J\f and the denominator D of eq. (4) has 
to simultaneously vanish. Hence by making M = T> = 0, 
we get the so called sonic point conditions (value of dy- 
namical and sound velocities in terms of horizon location 
r h ) as: 



Uh = a h 



1 



Ar h - 3 



(5) 



n 2 = f/£ 2 - \/£2 + £f, n± = (n 1 ±n 2 ) (8) 



(2) so that the three roots for rh comes out to be: 



X r h = - y + V+, 2 r h = - y - \ (n+ - iV3fL 



Th 



3 2V 



(9) 



However, note that not all V/j{i = 1,2,3} would be real 
for all "P 2 . It is easy to show that if \& > 0, only one root 
is real; if \& = 0, all roots are real and at least two of them 
are identical; and if 4' < 0, all roots are real and distinct. 
Selection of the real physical (r/, has to be greater than 
unity) roots requires the following discussions. 
Although we analytically calculate rh and other variables 
at rh, there is no way that one can analytically calculate 
flow variable at any arbitrary r. One needs to integrate 
eq. (6) only numerically to obtain the transonic profile of 
the flow for all range of r, starting from infinity and end- 
ing on to the actual event horizon. To do so, one requires 
to set the appropriate limits on Vi to model the realis- 
tic situations encountered in astrophysics. As £ is scaled 
in terms of the rest mass energy and includes the rest 
mass energy (£ = 1 corresponds to rest mass energy), 
£ < 1 corresponds to the negative energy accretion state 
where radiative extraction of rest mass energy from the 
fluid is required. For such extraction to be made possible, 
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accreting fluid has to posses viscosity and other dissipa- 
tive mechanism, which may violate the Lorenzian invari- 
ance. On the other hand, although almost any £ > 1 
is mathematically allowed, large values of £ represents 
flows starting from infinity with extremely high thermal 
energy, and £ > 2 accretion represents enormously hot 
(thermal energy greater than the rest mass energy) ini- 
tial (r — ► oo ) flow configurations, which is not properly 
conceivable in realistic astrophysical situations. Hence 
we set 1 £ ;$ 2. Now, 7 = 1 corresponds to isothermal 
accretion where accreting fluid remains optically thin. 
This is the physical lower limit for 7. 7 < 1 is not real- 
istic in accretion astrophysics. On the other hand, 7 > 2 
is possible only for superdense matter with substantially 
large magnetic field (which requires the accreting mate- 
rial to be governed by GR magneto-hydrodynamic equa- 
tions, dealing with which is beyond the scope of this pa- 
per) and direction dependent anisotropic pressure. We 
thus set 1 £ 7 £ 2 as well, so V2 has the boundaries 
1 {£,7} 2. However, one should note that the most 
preferred values of 7 for realistic BH accretion ranges 
from 4/3 to 5/3 0. 

Coming back to our solution for r^, we find that for our 
preferred range of V2, one always obtains \& < only. 
Hence the roots are always real and three real unequal 
roots can be computed as: 

\ h = 2V^cos (|) = 2v^cos 



where 



cos 



6 + 47T 



We find that for all 1 ^ V2 ^ 2, \h becomes negative. 
We then find that {V^, 3 ?^} > 1 for most values of our 
astrophysically tuned V2. However, it is also found that 
\h does not allow steady physical flows to pass through 
it, either u, or a, or both, becomes superluminal before 
the flow reaches the actual event horizon, or the Mach 
number profile shows intrinsic fluctuations for r < 77,. 
We obtain this information by numerically integrating 
the complete flow profile passing through T^. Hence it 
turns out that we need to concentrate only on V/, for 
realistic astrophysical BH accretion. Both large £ and 
large 7 enhances the thermal energy of the flow so that 
the accreting fluid acquires large radial velocity to over- 
come a only at the close vicinity of the BH. Hence rh 
anti-correlates with T>i. To obtain (dvjdr) and (da/dr) 
on AEH, we apply L' Hospital's rule on eq. (4) and we 
have: 



dr , 
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$ 12 = -$ 2 /2$i, $123 = V*2 - 4$i$ 3 /2$i 



$1 = 



Gr h (r h - 1) 
V4r h - 3 : 



$ 2 = 



4r ft - 3 



[4r h (7-I) -(3 7 -2)] 



8(^-1) 
(4r h -3) 5 



H (7 - l) 2 - r h (IO7 2 - 19 7 + 9) 



+ (6 7 2 -ll7 + 3)], $4 



2 {2r h - 1) - 7 (±h - 3) 
4 (r h - 1) 



(12) 

GR description of the sonic geometry for generalized 
acoustic metric is provided by Bilic 1]. Keeping in the 
mind that the surface gravity can be obtained by com- 
puting the norm of the gradient of the norm of the Killing 
fields evaluated at rh , we use the methodology developed 
in Bilic [l| to calculate the general relativistic expression 
for surface gravity, and Tah for our flow geometry is ob- 
tained as: 



Tah — 



fie 3 



rl (r h - 0.75) 



(a — u) 



4nk b GM BH [ (r h -l)-- j 

-jvhere the values of r^, {du/dr)h and (da/dr)h is obtained 
using the system of unit and scaling used in our pa- 
per. Stationarity approximation allows us not to dis- 
criminate between acoustic apparent and event horizon. 
Note again, that, since and other quantities appearing 
in eq (13) is analytically calculated as a function of V2, 
we provide an exact analytical value of the general rel- 
ativistic analogous Hawking temperature for all possible 
solutions of an spherically accreting astrophysical black 
hole system, which has never been done in the literature. 
If \JWh - 3(1/2 - l/$ 4 )($i2 - $123) > 1, one always 
obtains (da/dr < du/dr)h, which indicates the presence 
of the acoustic white holes at r^. This inequality holds 
good for certain astrophysically relevant range of V2, see 
following discussions. For a particular value of V2, we 
now define the quantity r to be the ratio of Tah an d 
T H , i.e., r = Tah/Th- t comes out to be independent of 
the mass of the black hole, which enables us to compare 
the properties of two kind of horizons (actual and acous- 
tic) for accreting BH with any mass. Our calculation is 
thus applicable for the widest possible range of BH can- 
didates, starting from the primordial holes to the super 
massive black holes at galactic cecntres. 

Figure 1 shows the dependence of t on V2 ■ The 
dependence is highly nonlinear. We plot r corespond- 
ing to the black hole solutions only (solutions for which 
(da/dr)h > (du/dr)h). Note that Tah > Th for some 
V-z . In figure 2., we scan the complete V2 space by de- 
picting its four important non-overlapping regions. For 
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Fig. 1: Variation of ratio of the analogous to the actual 
Hawking temperature (r) on conserved specific energy (£) 
and polytropic indices (7) of the flow. 
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Fig. 2: Parameter space classification for V2 for four 
different sets of values of Tah- 



accretion with Vi taken from the lightly shaded regions 
marked by < Tah < Th, the Hawking radiation dom- 
inates over its analogous counterpart. For high 7 high 
£ flow ('hot' accretion), AHR becomes the more domi- 
nant process compared to the actual Hawking radiation 
(Tah > Th)- For low £ and intermediate/high 7, acous- 
tic white holes appear (white region marked by WH, 
where da/dr < du/dr at 77,). For Vn belonging to this re- 
gion, one obtains outflow (outgoing solutions with Mach 
number increasing with increase of r) only. This has 
also been verified by obtaining the complete flow profile 
by integrating eq. (6). The dark shaded region (lower 
right corner of the figure) represents V2 for which is 
comes out to be real physical (r/i > 1) but $123 becomes 
negative, hence (du/dr)h and (da/dr)h are not real and 



Tah becomes imaginary. Note that both Tah > Th and 
white hole regions are obtained even for 4/3 < 7 < 5/3, 
which are the values of 7 for most realistic flow of matter 
around astrophysical black holes. Hence we show that the 
domination of the analogous Hawking radiation over the 
actual Hawking radiation and the emergence of analogous 
white holes are real astrophysical phenomena. 
However, the exact astrophysical nature of the acous- 
tic white hole solutions are not quite clear to us at this 
moment. We did not perform formal stability analysis 
on these solutions and it is possible that the white hole 
solutions are perhaps unstable, according to the recent 
proposal that the sonic white holes posses intrinsic insta- 
bility However, such detail investigations are beyond 
the scope of this paper and we would like to defer it 
for our future work. Note that calculations presented in 
this letter may be extended for astrophysical BH accre- 
tion discs as well. Bilic studied the relativistic acous- 
tic geometry for 2D axisymmetric space time and one 
can use such formalism to study AHR for GR and post- 
Newtonian multi-transonic, rotating, advective BH ac- 
cretion discs in the similar way as is done in this letter. 
Such calculations are in progress and will be presented 
elsewhere. 
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